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Abstract 

The production of charged hadrons and Kg mesons in the collisions of quasi-real photons has 
been measured using the OPAL detector at LEP. The data were taken at e'''e“ centre-of-mass 
energies of 161 and 172 GeV. The differential cross-sections as a function of the transverse 
momentum and the pseudorapidity of the charged hadrons and Kg mesons have been compared 
to the leading order Monte Carlo simulations of PHOJET and PYTHIA and to perturbative 
next-to-leading order (NLO) QCD calculations. The distributions have been measured in the 
range 10 < JE < 125 GeV of the hadronic invariant mass W. By comparing the transverse 
momentum distribution of charged hadrons measured in 77 interactions with 7 -proton and 
meson-proton data we hnd evidence for hard photon interactions in addition to the purely 
hadronic photon interactions. 
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1 Introduction 


Inclusive hadron production in collisions of quasi-real photons can be used to study the structure 
of photon interactions complementing similar studies of jet production in 77 collisions ||l|]. The 
photons are radiated by the beam electrons^ carrying only small negative squared four-momenta 
Q'^. They can therefore be considered to be quasi-real ~ 0) if the electrons are scattered 
at very small angles where they are not detected. For the “anti-tagged” event sample, events 
are rejected if one or both scattered electrons have been detected. 

The interactions of the photons can be modelled by assuming that each photon can either 
interact directly or appear resolved through its fluctuations into hadronic states. In leading 
order Quantum Ghromodynamics (QGD) this model leads to three different event classes for the 
77 interactions: direct, single-resolved and double-resolved. In resolved events partons (quarks 
or gluons) from the hadronic fluctuation of the photon take part in the hard interaction. The 
probability to hnd a parton in the photon carrying a certain momentum fraction of the photon 
is parametrised by parton density functions. 

We measure differential production cross-sections as a function of the transverse momentum 
and the pseudorapidity of charged hadrons and neutral Kg mesons. Since the distributions are 
fully corrected for losses due to event and track selection cuts, the acceptance and the resolution 
of the detector, they are directly comparable to leading order Monte Garlo models and to next- 
to-leading order (NLO) perturbative QGD calculations by Binnewies, Kniehl and Kramer p|. 
Until now, transverse momentum distributions of charged hadrons have only been measured 
for single-tagged events by TASSO [Q and MARK H 0] at an average (Q^) of 0.35 GeV^ and 
0.5 GeV^, respectively. We present the hrst measurement in anti-tagged collisions of quasi-real 
photons. Furthermore, the transverse momentum distributions in 77 interactions are expected 
to have a harder component than in photon-proton or meson-proton interactions due to the 
direct photon interactions. This will be demonstrated by comparing our data to the photo- 
and hadropro duct ion data measured by WA69 [^]. 

^Positrons are also referred to as electrons 
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At large transverse momenta (after crossing the charm threshold) the production of Kg 
mesons in photon-photon collisions is sensitive to the direct production of primary charm quarks 
in addition to the production of primary strange quarks, since the photon couples to the quark 
charge. Kg production in anti-tagged 77 collisions has previously been measured by TOPAZ 
and in single-tagged events by MARK II |^. 

In this paper, charged hadron and Kg production are studied using the full data sample taken 
in 1996 at e'''e“ centre-of-mass energies of 161 and 172 GeV corresponding to an integrated 
luminosity of about 20 pb“^. 


2 The OPAL detector 


A detailed description of the OPAL detector can be found in Ref. |0] , and therefore only a brief 
account of the main features relevant to the present analysis will be given here. 

The central tracking system is located inside a solenoidal magnet which provides a uniform 
axial magnetic held of 0.435 T along the beam axis0. The detection efficiency for charged 
particles is close to 100 % within the polar angle range | cos^l < 0.92. The magnet is surrounded 
in the barrel region (| cos^l < 0.82) by a lead glass electromagnetic calorimeter (ECAL) and a 
hadronic sampling calorimeter (HCAL). Outside the HCAL, the detector is surrounded by muon 
chambers. There are similar layers of detectors in the endcaps (0.81 < | cos6^1 < 0.98). The 
small angle region from 47 to 140 mrad around the beam pipe on both sides of the interaction 
point is covered by the forward calorimeters (FD) and the region from 25 to 59 mrad by the 
silicon tungsten luminometers (SW). From 1996 onwards, including the data presented in this 
paper, the lower boundary of the acceptance has been increased to 33 mrad following the 
installation of a low angle shield to protect the central detector against possible synchrotron 
radiation. 

Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector, a vertex drift chamber, a large volume jet chamber with 159 layers 
of axial anode wires and a set of ^ chambers measuring the track coordinates along the beam 
direction. The transverse momenta px of tracks are measured with a precision parametrised by 
cTpt/pt = y^O. 022 -|- (0.0015 ■ px)^ (pt in GeV/c) in the central region. In this paper “transverse” 
is always dehned with respect to the 2 ; axis. The jet chamber also provides measurements of 
the energy loss, dF/dx, which are used for particle identihcation 0]. 

The barrel and endcap sections of the EGAL are both constructed from lead glass blocks with 
a depth of 24.6 radiation lengths in the barrel region and more than 22 radiation lengths in the 
endcaps. The FD consist of cylindrical lead-scintillator calorimeters with a depth of 24 radiation 
lengths divided azimuthally into 16 segments. The electromagnetic energy resolution is about 
where E is in GeV. The SW detectors [|] consist of 19 layers of silicon detectors 
and 18 layers of tungsten, corresponding to a total of 22 radiation lengths. Each silicon layer 
consists of 16 wedge shaped silicon detectors. The electromagnetic energy resolution is about 
25%//F (F in GeV). 

^In the OPAL coordinate system the z axis points in the direction of the e“ beam. The polar angle 0, the 
azimuthal angle (j) and the radius r denote the usual spherical coordinates. 
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3 Kinematics and Monte Carlo simulation 


The properties of the two interacting photons {i = 1,2) are described by their negative four- 
momentum transfers Q^. Each is related to the electron scattering angle 9[ relative to the 
beam direction by 

Qi = -{Pi - P'i? - 2EiE'(l - cos 00, (1) 

where pi and p[ are the four-momenta of the beam electrons and the scattered electrons, respec¬ 
tively, and Ei and E[ are their energies. Events with detected scattered electrons (single-tagged 
or double-tagged events) are excluded from the analysis. This anti-tagging condition is met 
when the scattering angle 9' of the electron is less than 33 mrad between the beam axis and the 
inner edge of the SW detector. It dehnes an effective upper limit, Q^ax, values of Qj 

for both photons. The hadronic hnal state is described by its invariant mass W. The spectrum 
of photons with an energy fraction y of the electron beam may be obtained by the Equivalent 
Photon Approximation (EPA) [Q: 


f'i/eiy) 


a 


'1 + (1 


— log 




2 

min 



( 2 ) 


with a being the electromagnetic coupling constant. The minimum kinematically allowed 
negative squared four-momentum transfer Q^in is determined by the electron mass me- 


0‘^ ■ = 
Vmm 


2 2 
m^y 

1 - 2 /' 


( 3 ) 


The Monte Carlo generators PYTHIA |T^ and PHOJET |]^ have been used to simulate 
quasi-real photon-photon interactions. More details about the event generation can be found 
in Ref. |^. All possible hard interactions relevant to photon-photon interactions are included. 
The fragmentation is handled by JETSET |M. PYTHIA uses the SaS-lD parametrisation 


for the parton densities of the photon and PHOJET uses the GRV parametrisation |^. An 


approximation is used for the processes with primary charm quarks, i.e. where the charm quark 
is produced in the hard interaction. These processes are simulated using the matrix elements 
for light quarks. Subsequently the charm quarks are put on the mass-shell. 


4 Event selection and background 


The production of charged hadrons and Kg mesons was studied using the data taken at e’''e“ 
centre-of-mass energies, A/s^g, of 161 and 172 GeV with an integrated luminosity of about 
9.9 pb“^ and 10.0 pb“^, respectively. Photon-photon events are selected with the following set 
of cuts: 


• The sum of all energy deposits in the EGAL and the HGAL has to be less than 45 GeV. 

• The visible invariant hadronic mass, Hecal, calculated from the position and the energy 
of the clusters measured in the EGAL, has to be greater than 3 GeV. 
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• The missing transverse energy of the event measured in the ECAL and the forward 
calorimeters has to be less than 5 GeV. 

• At least 3 tracks must have been found in the tracking chambers. A track is required 
to have a minimum transverse momentum of 120 MeV/c, at least 20 hits in the central 
jet chamber, and the innermost hit of the track must be within a radius of 60 cm with 
respect to the z axis. The distance of the point of closest approach to the origin in the 
Y(f) plane must be less than 30 cm in the 2 ; direction and less than 2 cm in the ref) plane. 
Tracks with a momentum error larger than the momentum itself are rejected if they have 
fewer than 80 hits. The number of measured hits in the jet chamber must be more than 
half of the number of possible hits. The number of possible hits is calculated from the 
polar angle cos 9 of the track, assuming that the track has no curvature. 

• To remove events with scattered electrons in the FD or SW, the total energy measured 
in the FD has to be less than 50 GeV and the total energy measured in the SW has to be 
less than 35 GeV. These cuts also reduce contamination from multihadronic events with 
their thrust axis close to the beam direction. 

• To reduce the background due to beam-gas and beam-wall interactions, |(^o)| must be 
smaller than 10 cm where {zq) is the error-weighted average of the track’s z coordinates 
at the point of closest approach to the origin in the ref) plane. Beam-wall events with a 
vertex in the beam-pipe are rejected by requiring the radial position of the primary vertex 
in the ref) plane to be less than 3 cm. 


After all cuts 56732 events remain. 


All relevant background processes apart from beam-gas and beam-wall events were studied 
using Monte Garlo generators. Multihadronic events (e'''e“ —>■ qq(7)) were simulated with 
PYTHIA 5.722 |^. KORALZ 4.02 [|^ was used to generate the process e+e“ —> and 

BHWIDE to generate the Bhabha process e’''e“ —e+e“( 7 ). Processes with four fermions in 
the hnal state, including W pair production, were simulated with grc4f |T^, EXGALIBUR ||17|| , 
VERMASEREN and FERMISV . All signal and background Monte Garlo samples were 


generated with a full simulation of the OPAL detector [^. They were analysed using the same 
reconstruction algorithms as for the data. The main background processes are multi-hadronic 
e+e“ annihilation events and e+e“ —>■ events. Other background processes are found 

to be negligible. The multihadronic background is mainly reduced by the cut on the sum of 
the energy measured by the HGAL and the EGAL and by the cut on the missing transverse 
energy. The background from all these processes after the selection cuts amounts to less than 
1 %. 


The cut on the energy in SW and FD rejects photon-photon events with electrons scattered 
at angles 6' larger than 33 mrad and with an energy greater than 35 GeV in the SW or greater 
than 50 GeV in the FD. From the Monte Garlo, the rate of events with 9' > 33 mrad and 
energies less than 50 GeV is estimated to be negligible. The effective anti-tagging condition is 
therefore 9' < 33 mrad. 
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Analysis 


5.1 Correction procedure 

The measured transverse momentum and pseudorapidity distributions of the charged hadrons 
and the Kg mesons have to be corrected for losses due to the event and track selection cuts, for 
the acceptance and for the resolution of the detector. This is done with Monte Carlo events 
which were generated with PYTHIA 5.722 and PHOJET 1.05c. The data are corrected by 
multiplying the experimental distribution, e.g. of the transverse momentum px, with correction 
factors which are calculated as the bin-by-bin ratio of the generated and the reconstructed 
Monte Carlo distributions: 

V / generated / cl(T 

\dpT 

V ‘iPT / reconstructed 

As a correction factor the mean value from PYTHIA and PHOJET is used. The distributions 
of the pseudorapidity r] = — lntan(6'/2) are corrected in the same way. This method only 
yields reliable results if the migration between bins due to the hnite resolution is small. The 
bins of the pt and |p| distributions have therefore been chosen to be signihcantly larger than 
the resolution expected from the Monte Carlo simulation. The average transverse momentum, 
(Pt), and the average pseudorapidity, (|?7|), in each bin is calculated directly from the data, 
since detector corrections are small compared to the statistical errors. 

The visible invariant mass, HAis; is determined from all tracks and calorimeter clusters, 
including the ED and the SW detectors. An algorithm is applied to avoid double-counting of 
particle momenta in the central tracking system and the calorimeters [|I|. All distributions are 
shown for 10 < IE < 125 GeV where W is the hadronic invariant mass corrected for detector 
effects. To minimize migration effects when using Eq. § for the detector correction, the bins 
in W must be larger than the experimental resolution and the average reconstructed hadronic 
invariant mass, (ITrec), should be approximately equal to the average generated hadronic invari¬ 
ant mass, (ITgen)- The average (HAis) and the resolution on HAis as a function of the generated 
hadronic invariant mass ITgen are therefore shown in Fig. |l|a, where the vertical bars show the 
standard deviation (resolution) in each bin. The average (ITgen) as a function of HAis is plotted 
in Fig. 0b, where the vertical bars give the error on the mean. This plot is used to determine 
a correction function so that (ITgen)/ITrec ~ 1- The value of HAis measured in the detector 
is on average signihcantly smaller than ITgen- The relation between ITgen and HAis shown in 
Fig. 0b is almost independent of the beam energy and the Monte Carlo generator used. A single 
polynomial is therefore used to calculate HAec from HAis- The polynomial is obtained from the 
ht shown in Fig. 0b. It is applied to the data and the Monte Carlo. 

The efficiency to reconstruct photon-photon events in the detector, estimated by the Monte 
Carlo, is greater than 20% for ITgen > 10 GeV and greater than 60% for ITgen > 50 GeV. 
The trigger efficiency is dehned as the ratio of the number of selected and triggered events to 
the number of selected events. It was studied using data samples which were obtained using 
nearly independent sets of triggers. On average the trigger efficiency for the lowest IT range, 
10 < IT < 30 GeV, is greater than 97% and it approaches 100% for larger values of IT. Only 


(4) 


measured 
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lower limits on the trigger efficiency can be determined with this method and therefore no 
correction factor is applied. 


5.2 Charged hadron production 

For the charged hadron analysis only particles with a proper lifetime r > 0.3 ns are used to 
dehne the primary charged hadronic multiplicity in the Monte Carlo. The primary charged 
hadrons originate either directly from the primary interaction or from the decay of particles 
with a lifetime r < 0.3 ns including A and Kg decay products. The track selection criteria are 
dehned as in Section ^ In order to avoid regions where the detector has little or no acceptance, 
all measurements of charged hadrons were restricted to the range \ri\ < 1.5 (| cos^l^ 0.9). In 
this range, the resolution on px is given by CTp^/p^ ~ 0.02 (see Section |^) and the resolution 
on p by ~ 0.02. For the px distribution in the range 10 < IF < 125 GeV the correction 
factors as defined in Eq. ^ decrease from about 1.7 for px > 120 MeV/c to about 1.1 — 1.4 for 
Px > 2 GeV/c. The correction factor of about 1.6 for the rj distribution is nearly constant for 
IpI < 1.5. The PHOJET and PYTHIA correction factors differ by about 3 — 10 %. 


5.3 Kc production 


The Kg mesons are reconstructed using the decay channel Kg —tt+tt which has a branching 
ratio of about 69% [^. The reconstruction procedure is similar to the procedure described 


in Ref. [^]. It has been optimised to increase the efficiency for Ending Kg mesons in photon- 
photon events. Tracks of opposite charge are paired together. In addition to other quality 
cuts the tracks must have a minimum transverse momentum of 120 MeV/c and at least 20 jet 
chamber hits. The intersection of the tracks in the rcj) plane is considered as a secondary vertex 
candidate if it satisfies the following criteria: 


the radial distance between primary vertex and the intersection point must be greater than 
0.5 cm and less than 150 cm. For events with at least 6 tracks the primary vertex is fitted 
and for events with less than 6 tracks the beam spot reconstructed from tracks collected 


from many consecutive e'''e events during a LEP fill is taken as primary vertex 


• the difference between the radial coordinate of the secondary vertex and the radial coor¬ 
dinate of the first jet-chamber hit associated with either of the two tracks has to be less 
than 10 cm; 


• the radial coordinate of the tracks at the point of closest approach to the primary vertex 
has to be greater than 0.2 cm; 


• the angle between the direction of Eight from primary to secondary vertex and the com¬ 
bined momentum vector of the two tracks at the intersection point has to be less than 5°. 


In addition, a fit was performed for track pairs passing all these cuts, constraining them to 
originate from a common vertex. A correction procedure was used to compensate for the 
















energy loss of the pions in the inactive material of the detector. All secondary vertices satisfying 
|M(7r+7r“) — 0.4977 GeV/c^| < 0.02 GeV/c^ are considered to be Kg decay vertices, where the 
mass M is calculated assuming that both tracks are pions. Finally, the residual background 
is reduced by requiring at least 20 dE/dx hits. The two tracks are identihed as pions if the 
dE/dx probability for the pion hypothesis, that is the probability that the specihc ionisation 
energy loss in the jet chamber {dE/dx) is compatible with that expected for a pion, exceeds 
5%. 


In Fig. ^ the tt+tt” invariant mass M is shown for all identihed secondary vertices in the 
selected events before and after applying the dE/dx cuts. After all cuts the reconstruction 
efficiency for Kg —>• tt+tt” decays is about 35.5% and the purity is about 95.5% for pt(Ks) > 
1 GeV/c, |r/(Kg)| < 1.5 and 10 < IF < 125 GeV. 


6 Systematic errors 


The following systematic errors, common to the charged hadron and Kg measurements, are 
taken into account: 


• The correction factors are obtained using PHOJET and PYTHIA, separately. The re¬ 
sulting distributions are averaged to get the hnal result. The differences between the two 
distributions are used to dehne the systematic error. 

• The lower limit on the trigger efficiency is taken into account by an additional systematic 
error of 3% on the cross-section in the range 10 < IF < 30 GeV. 

• Systematic errors due to the modelling of the detector resolution for the measurement of 
tracks were found to be negligibly small in comparison to the other errors. The systematic 
error due to the uncertainty in the energy scale of the electromagnetic calorimeter was 
estimated by varying the reconstructed EGAL energy in the Monte Garlo by ±5%. 

• The limited statistics of the Monte Garlo samples, especially at large transverse momenta 
Pt, is also included in the systematic error. 

• The systematic error of the luminosity measurement is negligible compared to the other 
systematic errors. 


The systematic error of the Monte Garlo modelling and of the EGAL energy scale and, 
for the low IF region, the error from the trigger efficiency contribute about equally to the 
total systematic error. In the Kg reconstruction additional systematic errors were studied by 
varying the parameters of the secondary vertex hnder and the dE/dx cuts. The full difference 
between the results is used to estimate the contribution to the total systematic error from 
the Kg reconstruction, the Monte Garlo model dependence and the EGAL energy scale. The 
systematic error affecting the Kg reconstruction and the error from comparing the PHOJET and 
PYTHIA correction factors are of similar magnitude. The total systematic error was obtained 
by adding all systematic errors in quadrature. The total systematic errors are highly correlated 
from bin to bin. 
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7 Results 


The differential inclusive cross-section dcr/dpT for charged hadrons in the region {rjl < 1.5 is 
shown in Fig. ^ for different corrected W ranges together with the statistical and systematical 
errors. The corrected cross-sections are given in Tables |l| and ^ 


The measured differential cross-sections are compared to NLO calculations by Binnewies, 
Kniehl and Kramer [Q. The cross-sections are calculated using the QCD partonic cross-sections 
to NLO for direct, single- and double-resolved processes. The hadronic cross-section is a con¬ 
volution of the Weizsacker-Williams effective photon distribution, the parton distribution func¬ 
tions and the fragmentation functions of Ref. which are obtained from a fit to e^e“ data 


from TPC and ALEPH. The NLO GRV parametrisation of the parton densities of the pho¬ 
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is used with = 131 MeV and rric = 1.5 GeV/c^. The renormalization and factor¬ 
ization scales in the calculation are set equal to ^pt with ^ = 1. The change in slope around 
Pt = 3 GeV/c in the NLO calculation is due to the charm threshold, below which the charm 
distribution in the resolved photon and the charm fragmentation functions are set to zero. 


The cross-section calculation was repeated for the kinematic conditions of the data presented 
here at an average e’''e“ centre-of-mass energy = 166.5 GeV and for scattering angles 
6' < 33 mrad. For the differential cross-section dcx/dpT a minimum px of 1 GeV/c is required 
to ensure the validity of the perturbative QGD calculation. For the same reason the differential 
cross-section (ia/d\rj\ is restricted to the region px > 1.5 GeV/c. The scale dependence of 
the NLO calculation was studied by setting ^ = 0.5 and 2. This leads to a variation of the 
cross-section of about 30% at px = 1 GeV/c and of about 10% for px > 5 GeV/c. The 
NLO calculations lie significantly below the data for IF < 30 GeV for da/dpx and d(j/d|p|. 
The agreement with the data improves in the higher IF bins. The NLO calculation is shown 
separately for double-resolved, single-resolved and direct interactions. At large px the direct 
interactions dominate. It should be noted that these classifications are scale dependent in NLO. 

The Px distribution for 10 < IF < 30 GeV is compared in Fig. || to px distributions in yp 
and hp (h= 7r,K) interactions measured by the experiment WA69 |^. The hp data are weighted 
by WA69 in such a way that they contain 60% vrp and 40% Kp data to match the expected 
mixture of non-strange and strange quarks in the photon beam of the yp data. The WA69 data 
is normalised to the yy data in the low px region at px ~ 200 MeV/c using the same factor for 
the hp and the yp data. The px distribution of WA69 has been measured in the Feynman-x 
range 0.0 < xp < 1.0. The hadronic invariant mass of the hp data is IF = 16 GeV and the 
average (IF) is of similar size for the yp data. In the yy Monte Garlo the average (IF) is about 
17 GeV in the range 10 < IF < 30 GeV, i.e. the average values of IF in the different data 
samples are approximately the same. Whereas only a small increase is observed in the yp data 
compared to the vrp and Kvr data at large px, there is a significant increase of the relative rate 
in the range px > 2 GeV/c for yy interactions due to the direct process. A clear deviation is 
seen at large px from the exponential fall-off expected for purely hadronic interactions. 

The differential cross-section d(T/d|p| is compared to the predictions of the Monte Garlo 
generators PHOJET 1.10 and PYTHIA 5.722 in Fig. taking into account the anti-tagging 
condition 6' < 33 mrad. In PHOJET the suppression of the total yy cross-section is 
parametrised using Generalised Vector Meson Dominance (GVMD) and a model for the change 
of soft hadron production and diffraction with increasing photon virtuality is also included. 
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The photon-photon mode of PYTHIA only simulates the interactions of real photons with 
= 0. The virtuality of the photons dehned by enters only through the equivalent photon 
approximation in the generation of the photon energy spectrum, but the electrons are scattered 
at zero angle. This model is not expected to be correct for larger values of Q^. We have 
therefore simulated events with < 1 GeV^ with the photon-photon mode of PYTHIA and 
events with > 1 GeV^ and 6' < 33 mrad with the electron-photon mode of PYTHIA. 


The differential cross-section dcr/d|r 7 | shown in Fig. ^ is nearly independent of \ri\ in the 
measured range. The \r]\ distribution is reasonably well described by PYTHIA and PHOJET 
for pt > 120 MeV/c, apart from the high W region where PHOJET appears to be below the 
data. For transverse momenta pt > 1-5 GeV/c (Fig. ^ both Monte Garlo models underestimate 
the data signihcantly. The same behaviour is observed for the NLO calculation at low W, but 
the agreement of the NLO calculation with the data improves in the high W bins. The corrected 
cross-sections are given in Tables ^ and 


The differential inclusive cross-sections da/dpT and dcr/d|p| have been measured for Kg 
mesons with pT(Kg) > 1 GeV/c and |? 7 (Kg)| < 1.5. The pt and rj dependent cross-sections 
are presented in the W range 10 < IF < 125 GeV (Fig. ^ and Tables |]-|^). In addition 
the Pt distribution is shown for two separate IF ranges (Fig. ^ and Table |^). The results are 
compared to PHOJET and PYTHIA. Based on the PYTHIA simulation using SaS-lD, about 
half of the Kg are expected to be produced from charm quarks at large px where direct processes 
are dominant. Both Monte Garlo models signihcantly underestimate the Kg production cross- 
section in the low px region where most Kg are expected to originate from primary strange 
quarks. The distributions are reasonably well described by the NLO calculations which use 
the Kg fragmentation function htted to MARK H [^] and ALEPH data in Ref. |^. The 
change in slope between px = 2 and px = 3 GeV/c in the NLO calculation is again due to the 
charm threshold. The variation of the calculated cross-section as a function of px for different 
choices of scales, ^ = 0.5 and 2, is largest around the charm threshold, about 30 — 40%, and 
10 — 20 % elsewhere. 


8 Conclusions 


We present measurements of differential cross-sections as a function of transverse momentum 
and pseudorapidity for charged hadrons and Kg mesons produced in photon-photon collisions 
at LEP. The data were taken at e’''e“ centre-of-mass energies of 161 and 172 GeV. 

The differential cross-section da/dpx for charged hadrons is compared to NLO calculations. 
In the range 10 < IF < 30 GeV more charged hadrons are found at large px than predicted. 
Good agreement between the NLO calculation and the data is found in the highest IF range, 
55 < IF < 125 GeV. The Monte Garlo models PYTHIA and PHOJET both underestimate 
the cross-section for tracks with px > 1.5 GeV and |p| < 1.5 in all IF ranges. The shape 
of the differential cross-section dcT/d|p| is well reproduced by the NLO calculations and the 
Monte Garlo models. A comparison of the px distributions of the 77 data to px distributions 
measured in 7 p and ( 7 r,K)p processes at similar invariant masses shows the relative increase 
of hard interactions in 77 processes due to the direct component. The transverse momentum 
and pseudorapidity distributions of the Kg mesons are reasonably well reproduced by the NLO 
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calculations, but they are significantly underestimated by the Monte Carlo models PHOJET 
and PYTHIA. 
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10 < ID < 30 GeV 

30 < ID < 55 GeV 

Pt [GeV/c] 

(Pt) [GeV/c] 

da/dpT [pb/GeV/c] 

(Pt) [GeV/c] 

da/dpT [pb/GeV/c] 

0.12-0.28 

0.20 

(3.11±0.01±0.24)xl04 

0.20 

(1.06±0.01±0.08)xl04 

0.28-0.44 

0.36 

(2.70±0.01±0.20)xl0^ 

0.35 

(8.80±0.06±0.58)x103 

0.44-0.60 

0.51 

(1.60±0.01±0.10)xl0^ 

0.51 

(5.31±0.04±0.32)x103 

0.60-0.80 

0.69 

(8.00±0.06±0.46)x103 

0.69 

(2.67±0.03±0.14)x103 

0.80-1.00 

0.89 

(3.38±0.04±0.17)x103 

0.89 

(1.26±0.02±0.06)xl0=^ 

1.00-1.20 

1.09 

(1.48±0.02±0.08)x103 

1.09 

(5.95±0.12±0.22)x102 

1.20-1.40 

1.29 

(6.64±0.16±0.40)x102 

1.29 

(2.92±0.09±0.12)x102 

1.40-1.60 

1.49 

(3.29±0.11±0.21)x102 

1.49 

(1.55±0.07±0.06)x102 

1.60-1.80 

1.69 

(1.75±0.08±0.11)x102 

1.69 

(9.01±0.49±0.55)xl0i 

1.80-2.00 

1.89 

(1.00±0.06±0.07)x102 

1.89 

(5.96±0.39±0.38)xl0i 

2.00-2.20 

2.10 

(6.04±0.48±0.37)xl0^ 

2.09 

(3.29±0.29±0.13)xl0i 

2.20-2.40 

2.30 

(4.18±0.39±0.27)xl0^ 

2.29 

(2.15±0.23±0.20)xl0i 

2.40-2.60 

2.50 

(2.06±0.30±0.08)xl0^ 

2.50 

(1.64±0.20±0.07)xl0i 

2.60-2.80 

2.68 

(2.04±0.31±0.07)xl0^ 

2.70 

(1.01±0.16±0.08)xl0i 

2.80-3.00 

2.90 

(1.12±0.25±0.05)xl0^ 

2.88 

(9.18±1.55±1.10) 

3.00-3.50 

3.22 

(1.09±0.16±0.11)xl0^ 

3.21 

(4.26±0.65±0.68) 

3.50-4.00 

3.71±0.01 

(6.03±1.35±0.21) 

3.74±0.01 

(3.58±0.62±0.32) 

4.00-5.00 

4.36T0.01 

(2.99±0.82±0.33) 

4.38±0.02 

(1.07±0.25±0.04) 

5.00-6.00 

5.55±0.03 

(1.40±0.64±0.07) 

5.50D0.04 

(7.12±2.38±0.44)xl0-^ 

6.00-8.00 

- 

- 

6.94±0.07 

(2.74±1.22±0.15)xl0-i 

8.00-15.0 

- 

~ 

9.51±0.32 

(l.ll±0.73±0.36)xl0-^ 


Table 1; Differential inclusive charged hadron production cross-sections da/dpT for \v\ < 1-5 
and in the W ranges 10 < hh < 30 GeV and 30 < hh < 55 GeV. The hrst error is statistical 
and the second is systematic. No value is given if the error on (px) is less than 0.01. 
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55 < hh < 125 GeV 

10 < IV < 125 GeV 

Pt [GeV/c] 

(Pt) [GeV/c] 

da/dpT [pb/GeV/c] 

(Pt) [GeV/c] 

da/dpT [pb/GeV/c] 

0.12-0.28 

0.20 

(5.85±0.05±0.58)xl0^ 

0.20 

(4.78±0.02±0.38)xl04 

0.28-0.44 

0.35 

(4.75±0.04±0.46)x103 

0.36 

(4.05±0.01±0.31)xl0^ 

0.44-0.60 

0.51 

(2.92±0.03±0.28)x103 

0.51 

(2.42±0.01±0.17)xl0^ 

0.60-0.80 

0.69 

(1.53±0.02±0.12)x103 

0.69 

(1.21±0.01±0.07)xl0^ 

0.80-1.00 

0.89 

(7.25±0.15±0.47)x102 

0.89 

(5.33±0.04±0.27)x103 

1.00-1.20 

1.09 

(3.89±0.10±0.20)x102 

1.09 

(2.45±0.03±0.12)x103 

1.20-1.40 

1.29 

(1.91±0.07±0.09)x102 

1.29 

(1.14±0.02±0.06)x103 

1.40-1.60 

1.49 

(1.04±0.06±0.04)x102 

1.49 

(5.87±0.13±0.31)x102 

1.60-1.80 

1.69 

(6.20±0.41±0.18)xl0^ 

1.69 

(3.28±0.10±0.19)x102 

1.80-2.00 

1.89 

(3.78±0.33±0.20)xl0^ 

1.89 

(2.00±0.08±0.11)x102 

2.00-2.20 

2.09 

(2.53±0.27±0.06)xl0^ 

2.09 

(1.20±0.06±0.06)x102 

2.20-2.40 

2.28 

(1.72±0.22±0.04)xl0^ 

2.29 

(8.10±0.49±0.57)xl0i 

2.40-2.60 

2.50 

(1.10±0.18±0.02)xl0^ 

2.50 

(5.00±0.38±0.20)xl0i 

2.60-2.80 

2.69 

(7.89±1.49±0.23) 

2.69 

(3.81±0.34±0.18)xl0^ 

2.80-3.00 

2.92 

(5.09±1.15±0.34) 

2.90 

(2.70±0.29±0.11)xl0^ 

3.00-3.50 

3.23 

(3.18±0.59±0.04) 

3.22 

(1.77±0.16±0.07)xl0^ 

3.50-4.00 

3.70±0.01 

(2.07±0.48±0.10) 

3.72 

(1.16±0.13±0.04)xl0^ 

4.00-5.00 

4.46±0.02 

(9.30±2.44±0.58)xl0-i 

4.40±0.01 

(4.32±0.61±0.19) 

5.00-6.00 

5.38±0.04 

(3.30±1.67±0.64)xl0-i 

5.48±0.01 

(1.95±0.45±0.09) 

6.00-8.00 

6.81T0.09 

(1.24±0.76±0.05)xl0-i 

6.92±0.03 

(5.97±2.01±0.40)xl0-i 

8.00-15.0 

9.91T0.34 

(2.64±1.76±0.23)x10-2 

9.69T0.18 

(1.10±0.46±0.07)xl0-i 


Table 2: Differential inclusive charged hadron production cross-sections da/dpT for \v\ < 1-5 
and in the W range 55 < hh < 125 GeV and for all hh (10 < hh < 125 GeV). The hrst error is 
the statistical error and the second error is the systematic error. No value is given if the error 
on (pt) is less than 0.01. 
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10 < IT < 30 GeV 

30 < IT < 55 GeV 

\v\ 

G 

V 

) 

da/dr] [nb] 

G 

7] 

) 

da/dr] [nb] 

0.00-0.30 

0.15 

9.91±0.06±0.70 

0.15 

3.27±0.03±0.20 

0.30-0.60 

0.45 

9.98±0.06±0.71 

0.45 

3.33±0.03±0.19 

0.60-0.90 

0.75 

10.05±0.06±0.71 

0.75 

3.41±0.03±0.19 

0.90-1.20 

1.05 

10.01±0.06±0.71 

1.05 

3.43±0.03±0.20 

1.20-1.50 

1.35 

9.54±0.06±0.67 

1.35 

3.33±0.03±0.19 


55 < IT < 125 GeV 

10 < IT < 125 GeV 

\v\ 

G 

7] 

) 

da/drj [nb] 

G 

7] 

) 

da/dr] [nb] 

0.00-0.30 

0.15 

1.80±0.02±0.16 

0.15 

14.97±0.06±1.12 

0.30-0.60 

0.45 

1.85±0.02±0.17 

0.45 

15.16±0.07±1.13 

0.60-0.90 

0.75 

1.88±0.02±0.17 

0.75 

15.33±0.07±1.14 

0.90-1.20 

1.05 

1.93±0.02±0.16 

1.05 

15.35±0.07±1.15 

1.20-1.50 

1.35 

1.90±0.02±0.16 

1.35 

14.75±0.06±1.08 


Table 3: Differential inclusive charged hadron production cross-sections d(j/d|? 7 | for pt > 
120 MeV/c and in the W ranges 10 < hh < 30 GeV, 30 < hh < 55 GeV, 55 < hh < 125 GeV 
and for all W (10 < W < 125 GeV). The first error is the statistical error and the second error 
is the systematic error. 



10 < IT < 30 GeV 

30 < IT < 55 GeV 

\v\ 

(^ 

1 ) 

da/dr] [nb] 

(^ 

1 ) 

da/dr] [nb] 

0.00-0.30 

0.15 

(8.26±0.47±0.37)xl0-^ 

0.14 

(4.51±0.27±0.20)xl0-^ 

0.30-0.60 

0.45 

(9.33±0.49±0.49)x10-2 

0.45 

(4.36±0.27±0.09)x10-2 

0.60-0.90 

0.75 

(7.72±0.45±0.24)x10-2 

0.76 

(4.66±0.28±0.44)x10-2 

0.90-1.20 

1.05 

(7.95±0.47±0.35)x10-2 

1.05 

(4.39±0.29±0.27)x10-2 

1.20-1.50 

1.34 

(8.17±0.47±0.31)x10-2 

1.35 

(4.56±0.30±0.19)x10-2 


55 < IT < 125 GeV 

10 < IT < 125 GeV 

1^1 


1 ) 

da/dr] [nb] 

(^ 

1 ) 

da/dr] [nb] 

0.00-0.30 

0.15 

(3.07±0.24±0.07)xl0-^ 

0.15 

(1.62±0.06±0.09)xl0-i 

0.30-0.60 

0.45 

(2.91±0.23±0.09)x10-2 

0.45 

(1.65±0.06±0.08)xl0-^ 

0.60-0.90 

0.75 

(3.21±0.25±0.07)x10-2 

0.75 

(1.59±0.06±0.09)xl0-^ 

0.90-1.20 

1.05 

(3.14±0.25±0.12)x10-2 

1.05 

(1.55±0.06±0.07)xl0-^ 

1.20-1.50 

1.35 

(3.35±0.26±0.07)x10-2 

1.35 

(1.60±0.06±0.07)xl0-^ 


Table 4: Differential inclusive charged hadron production cross-sections d(T/d|77| for px > 
1.5 GeV/c and in the W ranges 10 < hh < 30 GeV, 30 < hh < 55 GeV, 55 < hh < 125 GeV 
and for all W (10 < W < 125 GeV). The first error is the statistical error and the second error 
is the systematic error. 
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Pt [GeV/c] 

10 < IK < 125 GeV 

(Pt) 

GeV/c] 

da/dpr 

pb/GeV/c] 

1.0- 

-1.2 

1.09 

±0.01 

206.2 

± 

17.4 

± 

16.1 

1.2- 

-1.5 

1.33 

±0.01 

100.2 

± 

8.9 

± 

8.5 

1.5- 

-1.9 

1.66 

±0.01 

32.9 

± 

4.5 

± 

3.7 

1.9- 

-2.4 

2.11 

±0.02 

13.5 

± 

2.6 

± 

1.3 

2.4- 

-3.0 

2.65 

±0.04 

5.2 

± 

1.3 

± 

0.7 

3.0- 

-4.0 

3.37 

±0.08 

2.0 

± 

0.7 

± 

0.2 

4.0- 

-5.5 

4.56 

±0.15 

0.4 

± 

0.3 

± 

0.1 


Table 5: Differential inclusive Kg production cross-sections da/dpT for pT(Kg) > 1 GeV/c and 
|77(Kg)| < 1.5 in the IV range 10 < hh < 125 GeV. The first error is the statistical error and 
the second error is the systematic error. 


\v\ 

(Ihl) 

da/d\p\ [pb] 

0.0- 

-0.3 

0.15 

± 

0.01 

61.7 

± 

7.2 

± 

5.2 

0.3- 

-0.6 

0.45 

± 

0.01 

63.1 

± 

7.3 

± 

5.4 

0.6- 

-0.9 

0.76 

± 

0.01 

72.4 

± 

7.8 

± 

6.1 

0.9- 

-1.2 

1.05 

± 

0.01 

70.2 

± 

8.0 

± 

4.7 

1.2- 

-1.5 

1.34 

± 

0.01 

58.0 

± 

7.5 

± 

4.0 


Table 6: Differential inclusive Kg production cross-sections d(T/d|? 7 | for pT(Kg) > 1 GeV/c and 
|? 7 (Kg)| < 1.5 in the W range 10 < hh < 125 GeV. The hrst error is the statistical error and 
the second error is the systematic error. 


Pt [GeV/c] 

10 < IK < 35 GeV 

35 < IK < 125 GeV 

(pt) [GeV/c] 

da/dpT 

[pb/GeV/c] 

(pt) [GeV/c] 

da/dpT 

pb/GeV/c] 

1.0- 

-1.2 

1.08 

±0.01 

127.3 

± 

14.6 

± 

12.6 

1.10 

±0.01 

75.0 

± 

9.5 

± 

5.1 

1.2- 

-1.5 

1.33 

±0.01 

73.6 

± 

8.4 

± 

6.6 

1.33 

±0.01 

30.0 

± 

4.3 

± 

2.3 

1.5- 

-1.9 

1.66 

±0.02 

19.8 

± 

3.4 

± 

2.5 

1.67 

±0.02 

13.1 

± 

3.0 

± 

2.3 

1.9- 

-2.4 

2.11 

±0.03 

9.8 

± 

2.8 

± 

1.3 

2.11 

±0.03 

4.4 

± 

1.2 

± 

0.4 

2.4- 

-3.0 

2.72 

±0.04 

1.7 

± 

0.9 

± 

0.2 

2.62 

±0.05 

2.9 

± 

0.8 

± 

0.6 

3.0- 

-4.0 

3.37 

±0.15 

0.9 

± 

0.4 

± 

0.3 

3.38 

±0.09 

1.0 

± 

0.6 

± 

0.2 

4.0- 

-5.5 


- 



- 



4.56 

±0.15 

0.3 

± 

0.3 

± 

0.1 


Table 7: Differential inclusive Kg production cross-sections dcr/dpT for pT(Kg) > 1 GeV/c and 
|? 7 (Kg)| < 1.5 in the W ranges 10 < IK < 35 GeV and 35 < hh < 125 GeV. The hrst error is 
the statistical error and the second error is the systematic error. 
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Figure 1: The relation between the generated hadronic invariant mass hFgen and the visible 
hadronic invariant mass hFvis for PHOJET and PYTHIA Monte Carlo events. The vertical 
bars show the standard deviation in each bin in (a) and the error on the mean in (b). The 
polynomial £t shown in (b) determines the correction function for hEvis- 



Figure 2: Distribution of the 7r’'“7r“ invariant mass M for all identified secondary vertices in 
the selected events before (open histogram) and after (hatched histogram) applying the cut 
|M(7r'''7r“) — 0.4977 GeV/c^| < 0.02 GeV/c^ and the dE/dx cuts used to identify Kg —7r’''7r“ 
vertices. 
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Figure 3: Differential inclusive charged hadron production cross-sections da/dpT for |? 7 | < 1.5 
and in the W ranges (a) 10 < hF < 30 GeV; (b) 30 < hF < 55 GeV; (c) 55 < hF < 125 GeV 
and (d) for all W (10 < hF < 125 GeV) measured at = 161 and 172 GeV. The data 
are compared to NLO calculations for > 1 GeV/c by Binnewies et al. together with the 
separate contributions of double-resolved, single-resolved and direct 77 interactions. The inner 
error bar shows the statistical error and the outer error bar the statistical and systematic errors 
added in quadrature. 
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Figure 4: The px distribution measured in 77 interactions in the range 10 < hF < 30 GeV is 
compared to the pt distribution measured in 7 p and hp (h= 7 r,K) interactions in the experiment 
WA69 ip. The hp and 7 p data are normalised to the low px region at px ~ 200 MeV/c. The 
cross-section values given on the ordinate are therefore only valid for the OPAL data. 
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Figure 5: Differential inclusive charged hadron production cross-sections d(T/d|77| for pt > 
120 MeV/c and in the W ranges (a) 10 < hF < 30 GeV; (b) 30 < hF < 55 GeV; (c) 55 < 
hF < 125 GeV and (d) for all W (10 < hF < 125 GeV) measured at = 161 and 
172 GeV. The data are compared to the PHOJET and PYTHIA simulation. The Monte Garlo 
distributions are plotted as histograms using the same bin width as for the data. The statistical 
error is smaller than the symbol size. The error bars show the statistical and systematic errors 
added in quadrature. 
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Figure 6: Differential inclusive charged hadron production cross-sections d(j/d|? 7 | for pt > 
1.5 GeV/c and in the W ranges (a) 10 < hF < 30 GeV; (b) 30 < hF < 55 GeV; (c) 55 < hF < 
125 GeV and (d) for all W (10 < hF < 125 GeV) measured at = 161 and 172 GeV. 
The data are compared to the PHOJET and PYTHIA simulation and to NLO calculations. 
The Monte Garlo distributions are plotted as histograms using the same bin width as for the 
data. The inner error bar shows the statistical error and the outer error bar the statistical and 
systematic errors added in quadrature. 
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Figure 7: Differential inclusive Kg production cross-sections (a) da/dpT and (b) da/dipl for 
PT(Kg) > 1 GeV/c and |? 7 (Kg)| < 1.5 in the W range 10 < hF < 125 GeV. The data are 
compared to the PHOJET and PYTHIA simulation and to NLO calculations. The data were 
taken at = 161 and 172 GeV. The Monte Garlo distributions are plotted as histograms 
using the same bin width as for the data. The inner error bar shows the statistical error and 
the outer error bar the statistical and systematic errors added in quadrature. 






Figure 8: Differential inclusive Kg production cross-sections da/dpT for pT(Kg) > 1 GeV/c and 
|p(Kg)| < 1.5 in the W ranges (a) 10 < hF < 35 GeV and (b) 35 < hF < 125 GeV. The data are 
compared to the PHOJET and PYTHIA simulation and to NLO calculations. The data were 
taken at = 161 and 172 GeV. The Monte Garlo distributions are plotted as histograms 
using the same bin width as for the data. The inner error bar shows the statistical error and 
the outer error bar the statistical and systematic errors added in quadrature. 
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